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In order to treat loops in the Lorentz-violating QED model, we present a derivation of the QED 
axial anomaly that specifically highlights the infrared origin of the effect. This is done using di¬ 
mensional regularization while treating 75 as a spectator. This enables us to revisit aspects in the 
fermionic sector of Lorentz-violating QED which have analogous structure to the axial anomaly. 

In particular, it is shown that both the Chern-Simons and photon mass operators are not induced 
in the one loop effective action. At lowest order in the Lorentz-violating parameter, we can de¬ 
fine a non-relativistic potential that captures the effects of vacuum polarization. This leads to a 
Zeeman-like effect in vacuum which lifts the two-fold degeneracy of the hydrogen atom S-orbitals. 


1 . INTRODUCTION 

Since the advent of quantum mechanics, the role symmetries play in constructing fundamental theories has been 
paramount. Lorentz invariance is the spacetime symmetry underlying our most successful theory of nature; the 
Standard Model of particle physics (SM). Despite its great successes, the study of possible violations of the symmetry 
structure of the model is a desirable endeavor. In general, the search for symmetry violations is one possible probe of 
new physics beyond the SM 0 , 0 . In partiulcar, Lorentz invariant theories with extra dimensions must incorporate 
a viable mechanism by which the symmetry gets spontaneously broken at low energies ii- An extra motivation 
comes from the proposal of emergent symmetries [sl-Ill]. One can envisage that the UV completion of the SM does not 
possess any of its observed symmetries and that the symmetry structure manifests itself only at low energies in the 
effective theory. In any scenario the symmetry is not expected to be exact and possible violations should be present 
even at low energies. 

Although the specific nature of symmetry violation is important on a fundamental level, the phenomenological 
implications of such violations can be systematically studied within an effective field theory (EFT) framework. This 
approach has opened the door for the construction of the Standard Model Extension (SME) , which is an 

extension of the SM to include all possible Lorentz-violating renormalizable operators composed from the SM fields and 
respecting the SM gauge symmetries. The theory nevertheless respects the so-called observer Lorentz transformations, 
namely invariance under the conventional boosts and rotations of special relativity. The violation comes about by the 
existence of tensor fields associated with the spacetime vacuum structure. These ultimately break rotational and boost 
invariance in any specific fixed frame which leads observables to depend on the orientation/speed of the experimental 
set-up. These tensors were shown to arise naturally in string theory 0 |. A similar approach was used in [Tj| to study 
possible violations of the U{1)em gauge symmetry of the SM. 

On the other hand, some classical symmetries may not survive after a theory is quantized. The most famous 
examples being the axial anomaly in the SM and the trace anomaly in conformal field theories. Although the study of 
anomalies remains an active area of research, by now the physics of anomalies is covered in most textbooks on quantum 
field theory. In perturbation theory, the traditional derivations tend to focus on the ultraviolet behavior of Feynman 
diagrams. However, through the use of dispersive techniques Dolgov and Zakharov [T^ followed by Frishman et al. 
M emphasized the infrared origin of the axial anomaly^. In the context of the trace anomaly, the same physics was 
elucidated in 03 through the explicit construction of the low-energy effective action that encodes the long-distance 
fluctuations of massless particles. With this physics in mind, we rederive the axial anomaly of massless QED (2ll - [^ 
using the more familiar dimensional regularization (DR). To highlight the key physics, no reference is being made to 
the definition of 75 in D dimensions. The axial anomaly should not be affected by the usual subtleties associated 
with 75 since it genuinely emerges from the infrared. Our derivation is instructive and, as we explain below, proves 
essential to discuss the Lorentz-violating QED model. 

The derivation of the axial anomaly opens the door to re-investigate certain aspects in Lorentz-violating QED 


^ See also II8H20II for dispersive derivations of anomalies that highlights the same physics. 
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(LVQED) that share identical structure with the perturbative analysis of the axial anomaly. Those features of the 
model were the subject of an enormous body of work, see for instance (23 - 13^ which is by no means a complete list. 
Primarily, we focus on the possible induction of the Chern-Simons and photon-mass operators in the one-loop effective 
action. We argue that the ability of our method to capture the infrared structure of the axial anomaly provides a 
thorough test for the physical significance of the previously obtained results. The consistent implementation of DR 
in the presence of 75 enables us to reliably tackle a rather peculiar feature of the Lorentz-violating QED, namely the 
radiative breakdown of gauge invariance which was first realized in ( 3 ^ l3^ . 

In general, the study of possible symmetry violations in electromagnetic phenomena is arguably the best way to 
constrain the parameter space of the SME. One typically expects that tree level processes lead to the tightest bounds, 
but in some cases loop effects might also become important offering in addition interesting insights about the structure 
of the theory. In this work, our focus is on a special operator in the fermionic sector of LVQED [T^ 

= I^QED - ( 1 ) 

where tp is a four component Dirac spinor representing any charged lepton or quark and is a spacetime-independent 
four-vector parameterizing Lorentz violation. The new term is CPT-odd since is invariant under CPT [l3| . 
Moreover, the U{1)em gauge symmetry of electromagnetism is manifest. 

We will consider the one loop effective action (EA) including the vacuum polarization effects. It was first noticed 
in (a that the symmetry-breaking term could potentially induce the Chern-Simons (CS) electromagnetic operator 
in the EA^. The CS operator reads^ 


^cs[A] = 


( 2 ) 


where is a constant four-vector proportional to b^. Nevertheless, the authors pointed out that the result is 
regularization-dependent. At 0{b^^), the structure of the vacuum polarization tensor is in close analogy with the axial 
anomaly graphs. However, the result is manifestly gauge-invariant and thus it is not possible to fix the ambiguity 
in the induced coefficient by a symmetry requirement. It was later argued in that a non-perturbative treatment, 
which we review in the body of the paper, could fix this ambiguity thus assigning a preferred value for the induced 
coefficient. Another striking feature of this particular treatment is the violation of the Ward-Takahashi identity which 
takes place at 0{b^) It is needless to mention that the loss of transversality threatens the consistency of the 

theory and its renormalizability. We critically review the non-perturbative formulation to argue against a particular 
aspect of the construction. 

After developing our DR-based method for the axial anomaly, we apply it to the vacuum polarization in LVQED. 
We find that both operators do not arise at one loop, in particular, gauge invariance is manifest. Our approach 
substantially differs from the previous literature in that it shifts the attention from the regularization issues to 
focus more on the physical content of the results. In particular, it enables us to raise and answer the question if 
these anomalous operators arise from the calculable low-energy physics. In addition, the consistent treatment of 
dimensionally regularized loops that contain 75 is essential to discuss the alluded to aspects of LVQED. 

Finally, the vacuum polarization contains finite pieces proportional to the Lorentz-violating parameter. At leading 
order in b^ and below the electron mass, higher derivative Chern-Simons-like operators are induced in the EA. 
This class of corrections have not received much attention in the literature. We find that they lead to distinct 
phenomenology which will be discussed in the context of the hydrogen atom. In particular, the Lorentz-violating 
background introduces a Zeeman-like effect in vaccum which splits the hydrogen atom ground state as follows 




8a^|b| 

Qtt Oq nP 


(3) 


where ag is the Bohr radius and m is the electron mass. This correction can be used to set a bound on b possibly 
stronger than the current bounds available in the literature. 

The plan of the paper is the following. In sectionj^l we re-derive the QED axial anomaly as the primary example of 
our method and to glean the main message of our analysis. In section|31 we review the non-perturbative formulation of 
the theory and point out a common inconsistency in the original computations which precisely lead to the appearance 
of the CS operator and the violation of the Ward identity. We then use our method to show that both effects do not 


^ The authors of 138113^ studied the theory over curved backgrounds. 

® The classical properties of such operator are very interesting and were discussed thoroughly in the seminal work of 0- Notice also that 
the CS operator changes by a total derivative under a gauge transformation, i.e. the Lagrangian density is not gauge-invariant. 
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arise, in particular, gauge invariance is manifest. In section|3]we use the non-vanishing part of the vacuum polarization 
at 0{b^) to compute the correction to the Coulomb potential which leads to a Zeeman-like effect in vacuum. Finally 
in section [5] we conclude and summarize our results. 


2. THE AXIAL ANOMALY USING DIMENSIONAL REGULARIZATION 

The original derivation of the axial anomaly plU^ is carried in 4-dimensions without the explicit use of a regulator. 
If on the other hand an attempt is made to use DR, one naively worries about handling 75 ^. In this section, we 
show that the axial anomaly is correctly reproduced assuming complete ignorance about the properties of 75 in D- 
dimensions. Anomalies in field theory are insensitive to the UV [13|, i-e. the result is unchanged by the unknown 
heavy physics that operates at short-distance. In particular, one should find this echoed by the insensitivity of the 
axial anomaly to the specific definition of 75 when DR is employed. In general, a regularization prescription is a 
deformation of the original theory in order to control the ultraviolet behavior. The exact prescription of 75 is part of 
the regularization procedure and thus one should expect to find the axial anomaly free of any dependence on the the 
latter. Indeed this is the case as we show next. 

Let us start by reminding that the classical Lagrangian of massless QED has a global chiral symmetry 

. (4) 

The Noether current associated with this symmetry reads 

J 5 = (5) 

and is conserved upon using the equations of motion. In perturbation theory, we can calculate the vacuum matrix 
element of the axial current to one-loop order 

( 0 |Jr(x)| 0 )= I q = k+p . ( 6 ) 

The amplitude is given by two diagrams and the second diagram is obtained from the first by merely interchanging 



{p,v) and (fc,A). Hence, it suffices to compute 



d°l f ^ A / . \ 

(27r )^{I — kY + P + {I + pY + iO j 


(7) 


where the photons are on-shell = 0 and the curly brackets denote a trace operation. We integrate over the 

loop momentum via conventional means and contract with to find the anomalous contribution to the divergence 
of the axial current 


=A‘^\k,p) + B’^\k,p) 


( 8 ) 


For a thorough review, see |4lL[SSl . 
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where 


le 


2 pi p1—x 


A'^^(k,p) I dx I dy 


/o ^0 


-In 


{x^ + (1 - y)f 


and 


■ (1 - - x|^)7‘'7a} 

^u\n ^_ {(^-iB{{x-l)^-yf)-i^{xl/i-yf)Y[x]/i + {^-y)f)} 

^ Jo ^ A 


and we defined 


A 2 11 

A = —q xy, - = - 7 + IndTT . 

e e 


(9) 

( 10 ) 


( 11 ) 


Notice that the tensor B'^^{k,p) is UV finite, and naively appear to have an infrared singularity. Nevertheless, this 
tensor is infrared safe as we show below. We now come to a central point of our derivation, namely that all UV 
divergences cancel without any reference to the definition of 75 in £)-dimensions. In other words, 75 is treated as a 
spectator. Let us isolate the divergent piece in the tensor and employ the usual identities 

7 “ 7 /^ 7 „ = (-2 + 2 e ) 7 ^, ( 12 ) 


along with the cyclic property of traces to find 

Adt. = [4/{7''^75} - 4A;''{7^^75} - 2p''{7^^75} + 2fc^{7''^75} + 2g''^{llif)o - fhb}] ■ (13) 

The above identities and trace cyclicity are true statements in Z?-dimensions. Notice first of all that traces with 
four Dirac matrices cancel out. Now the interesting feature is that the divergent tensor is anti-symmetric under the 
interchange (p, v) and (fc. A) and thus upon adding the two diagrams we find 

A':di,{k,p) + A^diApA) = Q ■ (14) 

At this stage the result is proven to be unambiguously finite and regularization-independent. In general, if one defines 
a certain extension of 75 then any trace becomes an analytic function of e 



where the tensors appearing in the expansion will depend on the prescription used for 75 . If the symmetry structure 
of the divergent tensor in eq. (TT^ was different such that eq. m is not satisfied then the final answer, despite being 
finite, would explicitly depend on the specific 75 -prescription. 

Now we can freely enforce the following limits on the traces 

lim{ 7 ^ 7 '' 7 “ 7 ^ 75 } = limfy'^y'^ys} = 0, p, z/, a, /? e 0,1,2,3 . (16) 

£->■0 e-J-O 

These limits must exist regardless of any 75 -prescription used, in particular, they ensure that the result of any finite 
matrix-element in the theory computed with DR is identical to the 4D result. We show next that indeed the anomaly 
is recovered using the above limits. 

Turning our attention to the finite part of A‘'^, we see that it comes comes from the 0(e) pieces in eq. (TT^ and 
the In A piece in eq. ®. Hence 

lim [A''\k,p) + A^-^ip, k)] = ^e^-^Pp^kp . (17) 

Moving to B''^{k,p), we find the expected feature that once the limits in eq. (fTBl) are enforced the potentially singular 
terms vanish identically and one finds ® 

lim [B’'\k,p) + B^’^ip, k)] = ^e^-^^Pp^kp . (18) 


^ Notice that this is indeed Avell-defined since the trace operation formally precedes performing the Feynman integrals. 
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Putting everything together yields the well-known result 

p2 r flD iD^ 

= . ( 19 ) 

This is one of our main results. Contrary to the naive expectation, we showed that using DR even in the presence of 
75 is innocuous if one is interested in an infrared effect such as the axial anomaly. Our method should be applicable 
to any computation sharing similar structure with the axial anomaly example. 


3. THE CHERN-SIMONS TERM AND THE WARD-TAKAHASHI IDENTITY 


In this section, we apply the method of the previous section to Lorentz-violating QED given by the Lagrangian 
in eq. H]). In particular, we focus on the vacuum polarization contribution to the one loop effective action up to 
quadratic order in the Lorentz violating vector . Our concern is two-fold; the value of the Chern-Simons coefficient 
and whether the Ward identity is violated. Before we do so, let us sketch the main ideas presented in the literature 
and point out an inconsistency in the original derivations. 

The first attempt to compute the Chern-Simons coefficient appeared in [2^ which was followed by plenty of discus¬ 
sions regarding the (un)-ambiguity of the induced operator, see for example (25l - l^ . In [3, the authors introduced 
a novel non-perturbative formulation of the theory to argue for a preferred non-vanishing value for the induced co¬ 
efficient. Non-perturbative simply means that one uses the exact fermion propagator rather than treating 
perturbatively as a new vertex in the theory. For example, the vacuum polarization (VP) tensor becomes 




d'^l f ^ i ^ _ i _1 

(27r)4 /-TO - 1 ( 75 "’' f 


( 20 ) 


In the conventional perturbative treatment the Chern-Simons coefficient is finite but ambiguous [l^ [^. To linear 
order in b^, there are two diagrams shown in figure with each being linearly divergent but the final result is finite. 
Nevertheless, one is free to shift the momentum variable in each diagram by different amounts and, due to the linear 
divergence, this renders the final answer undetermined. The issue here is very similar to the axial anomaly graphs, 
but with one subtle distinction: the answer, being proportional to the Levi-Civita tensor, is gauge invariant. In 
conclusion, there is no symmetry to impose on the matrix element which leaves the induced coefficient ambiguous. 



FIG. 2: The 0(b) contribution to the vacuum polarization. The gray square denotes the vertex —ifiys. 

On the contrary, there is a single diagram in the non-perturbative treatment given by eq. (EHl)- Practically, one 
expands the propagators in powers of b^ to perform the computation and evidently the same integrals at 0 {h) arises 
similar to the perturbative treatment. The new feature is that the momentum routing between the two integrals 
becomes correlated which gives a preferred value for the induced coefficient. Subsequently and using different 

regulators, a variety of answers were displayed in the literature [2^ - 1^ . A shocking aspect of this formulation is the 
violation of the Ward identity at 0[b^) [H, [33,13. 

Our primary goal in this section is to highlight an inconsistency in the original computations, in particular, the 
misuse of symmetric integration in divergent integrals. To keep our presentation self-contained, we review the cal¬ 
culations presented in [3 13 iii appendix El These results hinged on the use of symmetric integration by which 
logarithmic divergences cancel out and one can evaluate all integrals in AD without an explicit use of a regulator. By 
examining a well-known integral, we show below that symmetric integration is not a valid operation in the presence 
of divergences. We then move to apply our method to these calculations. 


3.1. The problem with symmetric integration 


Let us consider the following integral in AD 





(21) 
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This integral is linearly divergent and appears to vanish under a naive shift of the integration variable. Nevertheless, 
a careful treatment yields an unambiguous finite answer. We start by defining 

po = ipA, lo = ik, lo = ih ■ (22) 

We then Wick rotate to go to Euclidean space. To turn the difference in the integrand into a surface term, the second 
term is Taylor expanded 


{I - pY 






(23) 


By Gauss’s theorem, the AD integral is turned into a surface integral over a 3-sphere at infinity and thus all terms 
with two or more derivatives vanish. We find 


Ii{p) 


ipj lim 
l—¥00 




lilj 


(24) 


Spherical symmetry allows us to replace IJj = 5ijP/A and the integral becomes elementary. Transforming back to 
Lorentzian space yields the well-known result 


4 (P) = —Pm 


(25) 


On the other hand, one can arrange for this integral to vanish if symmetric integration is used instead. To see how 
this comes about, let us substitute the expansion of eq. (1251) in eq. (1211) and discard the quadratic pieces that do not 
contribute. Instead of using Gauss’s theorem, one can identically write 


Ii = ip^ 


dHf 


li 


Alilj 

J 


(26) 


Both pieces in the integrand are logarithmically divergent. If one allows symmetric integration via hlj 1%/A then 
clearly the result vanishes, in contradiction to the finite result in eq. (1251) . Since the latter is undoubted, we conclude 
that symmetric integration is not a well-defined operation for divergent integrals. In other words, a symmetric 
integration prescription in a divergent integral amounts to setting finite surface terms to zero. 


3.2. A vanishing CS coefficient and restoration of tranversality 


Based on the previous discussion, an explicit regulator must be introduced to carry out the computation. As we 
showed in the previous section, the use of DR is well suited for such computations. If the anomalous operators under 
study are genuinely induced by long-distance fluctuations, our method will indeed capture these effects. We will find 
that both operators vanish at one loop, in particular, gauge invariance is restored. 

Ideally, one might hope to perform the integral in eq. (1201) exactly and then expand the final answer in powers of 
b^. The complicated form of the propagators renders this task cumbersome, although an attempt was made in (^ . 
We choose to expand the propagators directly 


/-m-^75 / - m t - 


+ -A — {-itlY ■ 


-f 


t-' 


{-ihb) 




t — m 

{-ihY 


I-' 


(27) 


Evidently, the use of DR ensures the result matches onto the exact computation in the appropriate limit. Thus one 
finds 


where 


*nr(p) = (S'^'^(p) + E'^^(-p)) 


(28) 



/- 


m 


t + i>- 


m 


■ {-ihY 





(29) 
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Notice that a shift in one of the integrands has been used to cast the result in this form, which is indeed allowed 
by virtue of DR. It is convenient to separate into a potentially divergent piece and a finite one 


. (30) 

lOTT^ 

Let us analyze the potentially divergent piece, it reads 

= j dx (^ - In [(-1 + e)(l - + (-1 + e)(l - x){'y^Yfh 5 } 

+ - e (31) 

where use has been made of eq. (TT^ and A = rn? —p^x{\ — x). Using only the cyclic property of traces and the Dirac 
algebra, we find that terms multiplying 1 /e can be reduced to 


^ div. 


{phb} 


(32) 


In particular, all traces with four Dirac matrices cancel out. The above tensor is antisymmetric under interchanging 
/i and V while sending p —>■ —p, and hence is indeed finite. The development here is quite analogous to the axial 
anomaly. The finite parts are rather straightforward 

^^in. = J dx ^x (2 - 3x) In A + x‘^{2x — x^ — 1 )^ + - 

2 

/ dxx{2-x)'^ 

Jo ^ 

where the limits has been taken according to eq. dni). To arrive at the above result, one has to be careful with terms 
linear in m. These terms comprise traces of an odd number of Dirac matrices with 75 but since the integrals are 
logarithmically divergent, their traces can not be immediately dropped. Nevertheless, a careful computation shows 
once again that all divergences cancel identically and hence the whole contribution vanishes when the trace limits are 
enforced. The Chern-Simons operator is linear in p, and so we take the limit = 0 to find 

= 0 . (35) 

p2_).o 


(33) 

(34) 


There is no induced Chern-Simons term at one loop®. Notice also that DR guarantees the result is continuous in 
the massless limit as one can easily check by setting to = 0 at the start^. We notice here that the result does not 
suffer from any ambiguities related to regularization and is fully consistent with gauge invariance. 

Now we move to the next term in the expansion; the piece quadratic in 6 ^ which is the source of the violation of 
the Ward identity. The breakdown of the Ward identity posits a serious challenge to the consistency of the theory 
and its renormalizability. We set p = 0 from the start since the momentum-dependent terms necessarily satisfy the 
Ward identity js^- The polarization operator in eq. (OUll expanded to second order in b reads 

(36) 

where 


= 

= 


d^k -I- m){-i^j 5 ){ll: -|- 771)7“"(^ + + m)} 

(27r)'° (fc^ — rrF)'^ 

d°k {7^(1^ -I- + to)} 

(27r)'® (fc^ — rrfiY 


(37) 


We decompose the various tensors into divergent and finite pieces 

= D>^'' + -f (38) 


® This null result has been obained in l2flll . 

^ This does not take place for instance in the non-perturbative treatment as discussed in the appendix. 








with 


f ) t ^- = {7'^7“7"7^ (-*1^75)7" (-*I^75)7"}^/3.p 

+ 'y^l°‘{-ih5)Y{-ih5h^ + v}Jap + to'‘{7'"7‘'(-*^75)(-*I^75)}>/o 

F^^'' = m^{7^7“7'^(-i^75)7^(-i|i75)+7^7“7'^(-i^75)7^(-i^75) 

+ {-ih5){-ih5)l^ + v}Jap + rn^{'^^Y{-ih5){-ih5^)}Jo (39) 


and we defined 


JaPap - 24(47^)2 ('?“/5^P<^ + T]oia llPp + VPi^Vap) 

r _ * 1 r _ * 1 

12(47r)2'^“^TO2’ "^0 6(47r)2m4 



(40) 


Notice here that we do not use the AD property 7 I = 1. Also, we can ignore all terms with odd powers of m since they 
involve a trace of 7 f; an d an odd number of Dirac matrices. Focusing solely on the finite tensors, we easily reproduce 
the result found in 


2 

ie^ Imr {2b>^b'' + r^^^''b^) . (41) 

We can then use the identities in eq. m and manipulate the traces in the divergent tensors employing the cyclic 
property to find that the terms proportional to 1/e cancel identically similar to the previous calculations. The pieces 
proportional to e then give a finite contribution which reads 

2 

ie^ lim + D>^‘' + = + 1 ^ (26''6^ + V^''b^) (42) 

exactly canceling the piece in eq. (1411) . It is gratifying to see that gauge invariance is unambiguously manifest with 
our method. 


4. A CHERN-SIMONS-TYPE TERM AND THE HYDROGEN ATOM 


In this section we study the contribution to Ft^'^ which arises if ^ 0. Taking 11/“^ into account, we study the 
effect of the vacuum polarization on the spectrum of the hydrogen atom. The one loop vacuum polarization up to 
linear order in reads 


= (43) 

where 

7>2 r'^ A 

^^QED = ^ - P^P'') dx x(l -x)\n— (44) 

is the usual QED vacuum polarization tensor renormalized in an on-shell scheme. We start by computing the matrix 
element describing electron-proton scattering assuming the proton is a Dirac fermion. The amplitude reads 

iMeP = [u{k', s')j>^u{k, s)] [u{p', r')Yu{jp, r)] . (45) 

The different momenta and spins are labeled in figure Also, G^i, is the dressed photon propagator 




-tripp 


-^r|pc 


■m' 




(46) 


The dots indicate the iteration of the vacuum polarization tensor and q = p — p' = k' — k is the momentum transfer. 
The appearance of the Levi-Civita tensor renders the evaluation of the above summation cumbersome and so we 
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FIG. 3: Electron-proton scattering. The grey blob denotes the vaccum polarization insertion. 


choose to work at fixed order and truncate the series at one loop. The one-photon exchange gives rise to the famous 
fine structure correction 


E 


( 1 ) ^ Eil 

2m 



An \ 

i + i/ 2 j 


(47) 


where j is the total angular momentum quantum number. Now moving to the corrections due to vacuum polarization, 
we find the spin-independent Uhleng potential which reads 


^Uhleng ) 




,5(3) (r) 


(48) 


This potential is computed by expanding the vacuum polarization tensor in powers of /m? and retaining the lowest 
order piece which is a good approximation since the hydrogen atom wave functions are almost constant over the 
electron Compton wavelength [^. We do not display the shift in the spectrum due to the Uhleng term as it is 
irrelevant to our discussion. 

What we are after is the similar corrections that arise from flO"^. We work with spinors normalized as ■u(p, s)u{p, s) = 
1 and take the proton mass to be infinite. The photon propagator receives the following correction 




(49) 


where 




dx 


{3x^ — 2x) In A -t- {x^ — 2a:: -I- 1) — -t- a;(2 — x) — -- 


At low momentum transfer, we expand the above function 




127r2 


(50) 


(51) 


To streamline the discussion, we take the background vector to be space like, = (0,b). Hence, the amplitude 
simplifies to 


M P = 


1 


127r2m2 q2 


u{k', s')j^u{k, s)] [u{p', r')^''u{p, r)] 


(52) 


where the momentum transfer is purely spatial in the limit of infinite proton mass. Using some properties of the Pauli 
matrices, we end up with 


MIp=M^ + Ms 


(53) 


where 


M, = 


127r2m2 q2 


O’ ■ qq • b — q^cr • b 

2m 


$5® 


(54) 
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and the spin-independent piece reads 


M^ = 


1 (b X p) • q 




OTT^m^ 2m 

In the Born approximation, the potential is given by the inverse Fourier transform of the matrix element 

(fq 




(55) 


(56) 


This leads to the following potential 


V^{r) 

Vs{r) 


dSTT^m^r^ 


b • (p X r) 


2471^7 


4777-5 


• r b • r — 


4777-3 


b • I 


-a ■ b(5^^^(r) 


(57) 

(58) 


Using lowest order perturbation theory, it is easy to obtain the shift in the spectrum of the hydrogen atom. It is 
ideal to only consider the correction to the S-orbitals since the spin-independent piece in eq. dszj vanishes in this 
case®. In addition, the extra spin-dependent pieces due to QED corrections are suppressed by the proton mass (45| . 
The unique effect of the correction in eq. (l58ll is to lift the two-fold degeneracy of any S-wave, for example, the energy 
split in the ground state reads® 




8a^|b| 
Qtt Oq m® 


(59) 


where ag is the Bohr radius. Hence, the Lorentz-violating correction mimics the Zeeman effect. Torsion pendulum 
experiments [1^, set the following bound on the Lorentz-violating parameter®®: |b| < 0(1) x 10“^^en. Although 
the energy splitting in eq. (I59II could be used to set a direct and independent bound on |b|, we are unfortunately not 
aware of any existing experiment that can proble such an effect even indirectly. This will be the subject of future 
work. 


5. CONCLUSION 

We used dimensional regularization to discuss a simple rather instructive derivation of the QED axial anomaly. 
In particular, we showed that the usual difficulties associated with 75 do not matter in this case. The underlying 
physics being the infrared origin of anomalies in field theory. We carried over the intuition to re-investigate aspects 
of LVQED which caused controversy in the literature. We found that the Chern-Simons coefficient unambiguously 
vanishes and that the Ward-Takahashi identity is preserved. The latter is very important to ensure the theory is 
unitary and renormalizable. We argued that our method offers a reliable test to the results previously obtained in 
the literature. 

On the phenomenological side, we showed that the Lorentz-violating effects on the vacuum polarization give rise to 
Zeeman-like effect where the vector plays the role of a uniform magnetic field. Unlike the usual QED corrections, 
this lifts the degeneracy of the S-orbitals in the hydrogen atom. This unique feature should provide a window to set 
stringent limits on the Lorentz-violating parameter, a task we found tedious due to the lack of experimental input. 
Nevertheless, we would pursure this particular aspect in the future. Another interesting direction is to explore the 
predictions of the theory on gravitational phenomena especially in regard to the early Universe. 
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® Notice that the expectation value of (1/r®) diverges for 1 = 0 states, nevertheless, it makes sense to ignore the effect of this term based 
on symmetry. 

® We take the axis of quantization to be along b for simplicity. 

See also l48ll for bounds using atomic parity violation experiments. 
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Appendix A: Non-perturbative formulation 


In this appendix, we give a quick derivation of the results obtained in [^. In particular, we highlight the role 
of symmetric integration in obtaining finite results. As previously mentioned, an exact computation (to all orders in 
b^) was attempted in (^ . Since we are only interested in the consistency of symmetric integration, we will present 
the easier computation along the lines of |24| . For the convenience of the reader we adopt a similar notation despite 
taking a slightly different approach. 

= (Ai) 

where 




(A2) 


Explicitly we have. 






{ 7 ^(/ + m)Y{t + ^ + to ) 7 “ 75 (/ + f + m)} 
(27r)^ — m? + iO){{l + pY — w? + iOY 

d*l {7^(/ + m)7‘"(/ + m)7“75(/+^ + m)} 
(27r)4 (P — _|_ j 0 ) 2 ((l + p )2 — 77^2 _|_ jO) 


(A3) 


Recall that the above integrals do not represent two different diagrams: only global momentum shifts are allowed in 
the non-perturbative formulation. The remarkable feature here is that linear divergences cancel out at the integrand 
level as we show next. To simplify, we set p^ = 0 from the start and evaluate the traces to find 

= + (A4) 


where the divergent tensor reads 


^ = -4* 


/ 


dH {P + 21 ■ + 2e^>^'^''PlaPfi 

(27r)"‘ {P — m'^ + iO){{l + pY — m'^ + iOY 


and 


_ _Y_ AiPva 

4n. -g^ 2 ^ PP 


(A5) 


(A6) 


We want to isolate the linear divergence prior to performing any loop shifts. We start by manipulating the divergent 
tensor 


Similarly, 


and 



(27r)4 [P — mP + *0)((Z -|- pY — rrP + zO)^ 

dH 

{2'kY {P — rrP + zO)((Z + pY — rrP + zO) 



dH + Pe'^^f^^pp + 

(27r)"‘ {P — rrP + z0)2((l -|- pY — nrP + zO) 

dH 

{27rY {P — m'^ + iO){{l + pY — + iO) 


jtiva _ tJ.Pua 

4n. -g^ 2 ^ PP 


(A7) 


(AS) 


(A9) 


Inspection of eqs. (IAtI) and (IA8I) shows that the linear divergences cancel out at the integrand level and the remaining 
divergences are at most logarithmic. The computation proceeds by introducing Feynman parameters and performing 
the loop integral. So far we have not introduced any kind of regulators, and here is where symmetric integration kicks 
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in. Using the latter prescription, one readily finds that the logarithmic divergences cancel out at the integrand level 
as well. After this cancellation, we can collect the finite pieces to find 

= (AlO) 

identical to the result in Notice that if the fermion mass was set to zero at the start, the result becomes 

= = . (All) 

This discontinuity alone is enough to suspect the validity of the procedure. A similar story takes place for 11^2 {Q) [3^ 
with the end result displayed in eq. m- We do not include the full derivation since it is very similar to the above 
computation. 
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